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I. INTRODUCTION 


Q 


There has been a lot of interest in understanding the spin-3/2 (or the Rarita-Schwinger 

This spin-3/2 field, or the 
his field in 
3| and j4| the 


) fields after the introduction of the supergravity theory 
gravitino, is the superpartner of the graviton. However, the consideration o 
curved spacetimes, particularly in black hole cases, is relatively sparse. In Refs, 
gravitino field was analyzed in the Kerr and Kerr-Newman backgrounds. The authors there 


used the Newman-Penrose formalism , following the earlier work of Chandrasekhar on the 
Dirac field jb]. The Newman-Penrose formalism is very useful in dealing with perturbations 
of the black hole in four dimensions, but the extension of the formalism to higher dimensions 
is not straightforward j^. For this reason we try, in this paper, to implement a novel and 
more direct approach to considering the gravitino fields in spherically symmetric spacetimes. 
In this approach we shall use the eigenspinors and the eigenspinor-vectors on spheres to 
separate out the angular parts of the gravitino fields. We can then obtain the radial equations 
for the various held components. 

These held components are gauge dependent in general. The gauge symmetry here is that 
of supersymmetry. In the Newman-Penrose formalism the radial equations for some scalar 
potentials are considered, where these potentials are gauge invariant. In order to compare 
with this formalism we also need to construct gauge invariant combinations out of our held 
components. The radial equations of these gauge invariant variables can then be used in 
the analysis of the evolution of spin-3/2 perturbations in spherically symmetry spacetimes, 
like the Schwarzschild background we concentrate on in this paper. We hope to extend this 
formalism to higher dimensions in our subsequent works. 

The study of the evolution of small perturbations in black hole backgrounds is an old and 
well established subject js, 9|. Actually, we know that this evolution, at intermediate times 
is dominated by damped single frequency oscillations. These characteristic oscillations have 
been termed quasinormal frequencies, and depend only on the parameters characterising the 
black hole. In this sense we can say that black holes have a characteristic sound 10|. 

Note that since the discovery of quasinormal oscillations of black holes, many studies have 
been done on quasinormal modes (QNMs) of various spin helds and a considerable variety of 


ana. 


11 


ytical, semi-analytical and numerical methods have been developed to determine them 


15| . However, an interesting problem is to consider the study of the QNMs of spin-3/2 
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fields, where the motivation for studying spin-3/2 QNMs from black holes in this paper is 
two-fold; the first of these being from the theoretical point of view, where the scarcity of 
work done in four dimensions Q, H] (and in a future work, greater than four dimensions 
[l^ ) is an omission in the literature. Our calculations serve to hll this gap. 

Secondly, having a complete catalog of all QNMs would be a necessary precursor to 
eventually studying the emission rates of collider produced, or TeV scale black holes |l9 |. 
Especially given that in super gravity theories the gravitational field is coupled to the mass¬ 
less spin-3/2 Rarita-Schwinger held, that acts as a source of torsion and curvature |20l. |21|. 


Furthermore, in many supergravity models, the supersymmetric breaking scale is related to 
the gravitino. As such, the gravitino could be the lightest supersymmetric particle, or even 
the next to lightest supersymmetric particle (which could live long enough to be directly 
detectable in collider environments). As such, spin-3/2 helds could have signihcant effects 
on phenomenology. 

Other than the quasinormal modes one is also interested in the greybody factors and 
the emission cross-sections of Hawking radiation of various spins when black hole evolution 
in collider environments are studied. The related absorption probabilities of the gravitino 
scatterings can be readily considered using the radial equation mentioned above. 

The approach we shall take is to use the WKB method, as it allows a more systematic 
study of QNMs, as well as the absorption probabilities, than outright numerical methods. 
We will use this, and previous results with WKB and similar methods , to compare with 


the newer method developed by some of us, whic 


the asymptotic iteration method (AIM) 


3 , 22 


1 can be more efficient in some cases, called 


24l |. In the case of absorption probabilities 


we shall also study the low energy case using the Unruh method. 

As such, this paper shall be structured as follows. In the next section we shall use 
the eigenspinors and the eigenspinor-vectors on the 2-sphere to separate out the angular 
parts of the gravitino, or the Rarit a-Schwinger, equations. In Section III, we construct 
the gauge-invariant variable of the gravitino field. The corresponding radial equations are 
compared with those obtained using the Newman-Penrose formalism. In Sections IV and V, 
the quasinormal mode frequencies and the absorption probabilities of the gravitino field are 
studied using the radial equation of the gauge invariant variable, respectively. Conclusions 
and discussions are presented in the last section. Some properties of the eigen-spinors and 
the eigen-spinor vectors relevant to our discussions are given in the Appendix. 
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II. MASSLESS RARITA-SCHWINGER EQUATION 


We start with the massless Rarita-Schwinger equation 

= 0 , ( 1 ) 

where we have used the notation 

_ ^ 2 ) 

Here we consider the Rarit a-Schwinger held in a Schwarzschild black hole background with 
the metric 


ds^ = —f{r)dt^ + /(r) ^dr"^ + r‘^dVt\ , 


( 3 ) 


where /(r) = 1 — 2M/r. Note that from now on an overbar represents quantities on the 
2 -sphere. 

Our choice of Dirac matrices are as follows: 


(g) 1 ^ y = —{ia^ (g) I) , 

v/ 

y = (T^ (8 y ^ = -((T^ o 7''^*) , 

r 

7 ^ = (^1 ^ Y = ® 1 ) ) 


( 4 ) 


where 1 is the 2x2 unit matrix, a* {i = 1, 2, 3) are the Pauli matrices, and 7 * {i = 1, 2) are 
the Dirac matrices for a 2-sphere. The corresponding non-zero spinor connection components 
are 

rt = -y<T' 0 i), 0 = 0 , 

o, =i0f„.+y(i<T=0 7«,). (5) 

Now we can work on Eq. ([T]). First, the product of Dirac gamma functions are found to 
be 


7 


tdiV 


7 


rdiOj 




( 6 ) 
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where Then we separate the helds into their r-t and angular parts. Since '0t 

and '0r behave like spinors on the 2-sphere, we can write 


V't = 0t (g) , 'tpr = (pr ® , 


(7) 


where ipx is a eigenspinor with eigenvalue i\. For completeness the eigenspinors on the 


2-sphere are briefly described in the Appendix. 0^ and 0^ are functions of r and t. Since 
the Schwarzschild metric is static, one can consider one frequency at a time and assume the 
time dependence to be 

00 . behaves like a spinor-vector on the 2-sphere. We can therefore write 


= 00^^ (8 V0,0A + «) 'yei^x , 


( 8 ) 


where we have represented the two sets of eigenspinor-vectors on the 2-sphere by 'Vg.ilJx and 


70i0A- Some properties of these eigenspinor-vectors are also presented in the Appendix. 

The equations of motion for 0^, 0^ and 00. can be derived by analyzing Eq. ([T]). For 
simplicity, we shall take the Weyl or temporal gauge such that 0t = 0. Gauge symmetry of 
the massless Rarita-Schwinger equation will be discussed in more detail in the next section. 
For /i = t in Eq. ([1]), we have 



(9) 


Using the definitions in Eqs. (I7j) and (IH]), Eq. ([9]) can be simplified to 




( 10 ) 


Hence, we have the first equation involving 0^, 0^^^ and (j)f\ 



y7(!Cr“)) <l>r 


0 . ( 11 ) 


Similarly for p = r we have from Eq. (ITI) , 
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Again using Eqs. ([7]) and ([8]), Eq. flT^ becomes 


■ 4/f 


ID. “AZffUi'i, 

2 r 'D» + 4,/7‘^ 'Do 2r 




'/7,U(2I + 

r 


6» T —1^ v^e T + “T*^ 'f'e 


-^ = 0 . 


(13) 


Note that we have taken the Weyl gauge: <pt = 0. Hence, we have the second equation of 
motion 


l\^/J ^ i\f ^ 1 2 ’\ Ab 

77 2r A^/J 2r 


\ff r 

Lastly, for n = Oi'm. Eq. ([1]), we have 


. ^\/7 f 

Ot H- (T + ^—F=(T H-cr 


2 1 J.(2) 


2/7 


= 0 . 


(14) 


= 0 


-a 0 ^ + 

r V’^V/ 


(i/)at + —a'^dr + -^j=aA 
r Ar^ ) 




. / 1 /' 1 ^ *A/7 2 

9t-7r^ - -rcx +-a , 

/r 4 r / \ TyJ j 


*A .. 3. „ zA\/7 20 *A/ 2 \ 1(1) 

-(*a )ai + ^—a a, + ) 0 / 




7^*0 = 0 . 


(15) 


Therefore we have two more equations of motion 


/7 


/7 


/' 


+ I — j={ia-^)dt H-+ 7 - 7 = 0 “ 

r V rv J ^ 4rv J 


2 1 ^(b 


= 0 , 


(16) 


and 


„ / 1 /' 2 


iXf 




21 Ab 


{i(j‘^)dt + ^^a'^dr + . ^ 0 / = 0 . 


^/7 


4r/7 


(17) 


In summary we have four equations of motion, Eqs. dm), dm, dm, and dm for three 
functions 0^, (j)^g\ and (j)f\ One can show that of these four equations, only three of them 
are independent. On the other hand, as we shall see below, the functions 0^, \ and 0^^^ 
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are not gauge-invariant. Therefore, in the next section we shall consider the gauge symmetry 
of the Rarita-Schwinger equation in detail. Then we shall construct the appropriate gauge- 
invariant variables and derive their corresponding radial equations of motion. 


III. GAUGE-INVARIANT VARIABLE 

In this section we consider the gauge freedom for Eq. ([T]). Consider the transformation 
of V'm, 

V’l, = , (18) 

where is a Dirac spinor. Plugging this into Eq. ([1]), we must have 

= 0 (19) 

for the equation to possess this gauge symmetry. Now, 

= ^7^"“^.ap.7"7> • (20) 

Note that to obtain the result for the commutator, [Vy, Vq], we have assumed that there 
are only gravitational interactions. For example, for a charged black hole background there 
should be extra terms proportional to the electromagnetic field strength. 

To simplify this expression in Eq. fl20|) further we need the following identity. From the 
symmetry of the Riemann tensor. 


+ Ruatifi + Raiiufi) = 0 ^ T R^iuaji + ‘^l°‘Rap) = 0 

^ 7^7"7“i?M-/3 = -27“^a/3 • (21) 

Also we have 

rrr-f^R,ua0 = -2R . ( 22 ) 

Using these identities Eq. fl20|) becomes 

= ^{2rR^^-YRW- (23) 
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This is zero for Ricci flat spacetimes like the four dimensional Schwarzschild spacetime that 
we look at below. However, even for de Sitter and anti-de Sitter spacetimes it does not 
vanish, so we need to modify the covariant derivative in those cases in order to respect the 
gauge symmetry. 

Consider the gauge transformations on <pt and 0^- Take (p = 0 $§ then Eq. flTS]) 
becomes 

= + Vtv? ^ (t)[ = (l)t + dt(j) - , (24) 

1pr + Vr'-P ^ 0), = • (25) 


The gauge transformation of the angular components of is more complicated. 
V'e, = iJOi + 

^ cl)f^ ® ® 7ei0A = + 0 ) <8 Ve,0A + ^06^^ + (8) 


^ = 0 ^^ + 0 ; 0 '^^^ = 0 ^ + ^(*^^)0 • 

We see that indeed 0,-, 0^, and .^7 tne not gauge invariant. 


(26) 


4 = + 


Considering the gauge transformations of these functions, one can construct gauge- 
invariant combinations. For example, a simple one would be 

(27) 

In the following we shall look for the radial equation for $ and then we shall show that it is 
the same as the one obtained in Refs. and jj]. 

Now we go back to the equations of motion in Eqs. ffTT]) . ffTTD . ffTB)) . and (ITT)) . Only three 
of them are independent and we shall consider Eqs. fllll) . (flT)) . and (II6p . From these three 
equations we can derive the equation of motion for d). 


^/f + Acr^ j 


+ + ^(v7-Aa3)<I) 




= {^/f-XcT^) /0.4> + |^(y/7-Aa3)$ 


(28) 


Note that in arriving at this form we have used the fact that / = 1 — 2M/r. 

Next we try to transform the radial equation, Eq. (|28ll . into the Schrodinger form. To do 
that we decompose the equation into its components, 

$2 


$ = 


(29) 


8 






and Eq. fl28|) becomes two equations 


fri) - V (^) • 

where we assume the time dependence of the helds is 

To further simplify the equations we make the following transformation, 


(30) 

(31) 


*1 = 




,\/7+A 

Then Eqs. flSO]) and (ITT]) become 
d 


<hi 


$5 — 


yj-A 


< 1 > 


2 • 


(32) 


dr* 


— IE $1 = 


dr* 


+ IE I $2 = ioj^i 


(33) 


with d/dr* = f d/dr and 


r ^ 2 -/ 


(i “ 5 ) U + 5 ) (i +1) v7 


(i+r-/ 


(34) 


where for convenience we have written |A| = j + 1/2 with j = 3/2, 5/^. .. 

Eq. (ITT]) is in the form of the so-called SUSY quantum mechanics 2^ with IE being the 
superpotential. Now we can construct two Schrodinger-like equations. 


d^ ~ 

—$1 + Ei$i = 

dr; 


d 2 


dr* 


$2 V 2$2 — *^2 


(36) 


where the isospectral potentials are 


E,2 = ±/^ E . 
dr 


(36) 


More explicitly. 


Ul ,2 = 


(a - 1) (a + 1) (a +1 ) \AT 
'■M(A-i) (a + I) + E]' 


2 M 


X 


±- 


2M^ 


J - 


J + 


J + 


^ 2M ^ 3M 


Tl 


(37) 


El is the same potential obtained in 3|, 1^. To see the behavior of Ei for various values of j 
we have plotted them in Fig. [D for 7 =3/2 to 11/2. They are similar to the barrier potentials 


of perturbations of other spins 


10 


26|. 
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FIG. 1: Effective potentials of the gravitino field with j = 3/2 to j = 11/2. 

IV. QUASINORMAL MODES 

In this section we consider the QNMs corresponding to the potential Vi for the gravitino 
in a Schwarzschild spacetime. For notational simplicity, we shall write V instead of Vi. 
The calculation of the QNM frequencies was pioneered by Chandrasekhar and Detweiler 
ll| using numerical methods. Due to the particular boundary conditions for a QNM, the 
numerical code is not easy to implement. In order to obtain the QNM frequencies in a more 
efficient manner, other semi-analytical methods have been devised, notably the continued 
fraction 12 1 and WKB methods jl3. 141. 

Up until recently, the WKB method has perhaps been the most popular method to 
evaluate the black hole QNM frequencies. The method was put forth by Schutz and Will 
3, then extended to third order by Iyer and Will [3 > and even to sixth order by Konoplya 
271]. In fact, the lowest order in the WKB approximation can be viewed as the large angular 
, that is, the large j limit in the potential of Eq. fl3711 . In this limit, the 


momentum limit 


QNM frequencies are given by the formula 




(38) 


where n = 0,1, 2,... is the mode number and Vq is the maximum of the potential. For 
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j —oo, the potential becomes 


^1 


j^oo 




(39) 


Note that from here on we shall take M = 1. The maximum of the potential in this limit is 
located at r = 3. Hence, 

2j2 


Hn 


01 j—>oo 


L. 

27 


V"\ ■ 

^0 Ij-i-oo 


729 


(40) 


The QNM frequencies in this limit are therefore 

1 


00 


j-^co 


J - iin + 2 ) 


(41) 


3\/3 L 

To obtain the QNM frequencies in small values of j, especially the fundamental ones, it 
is necessary to go higher orders in the WKB approximation. Using the formula put forth 


in Ref. 


14| for third order and in Ref. 27| for sixth order of the WKB approximation, we 


have evaluated the frequencies for j = 3/2 to 13/2, and with n < I {I = j — 3/2) since the 
WKB approximation is known to be accurate for low-lying modes. These values are listed 
in Table Hand can be compared with the ones given recently in Ref. 16|. 


In a series of papers 


[ 15,22 


24| we have developed a new semi-analytic method to evaluate 


the QNM frequencies called the asymptotic iteration method (AIM) as an alternative to the 
WKB approach. To implement the AIM approach in Eq. (1351) . we hrst make a coordinate 
transformation y'^ = 1 — 1/r, so that the region 1 < r < 00 is now in a compact one 
0 < 1 / < 1. Then we extract the asymptotic behaviors of <hi and write 


i>, = '{/-“(l - ■ 


(42) 


The function xiv) satishes the equation 


x"(y) = + sa(y), 


(43) 


where 

8ioo{l-Ay‘^ + 2y^) I - 5y^ 

° - y{l-y^y y{l - y^) ’ 

64a;^(2 — y"^) 32ioo 

(1 _ y2y 1 _ ^2 

(2j - l)(2j + l)(2j + 3)[(2j + 1) - 2yf[l + (2j + l)y + y^] 
2/(1 - 2/^)^[(l - - 4j2) 4|/2 - 2y^f 


(44) 


(45) 
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TABLE I: Low-lying (n < I, with I = j — 3/2j gravitino quasinormal mode frequencies using the 
WKB and AIM methods. 




WKB 

AIM 

1 

n 

3rd Order 

6th Order 

150 iterations 

0 

0 

0.3087 - 0.0902i 

0.3113 - 0.0902i 

0.3108 - 0.0899i 

1 

0 

0.5295 - 0.0938i 

0.5300 - 0.0938i 

0.5301 - 0.0937i 

1 

1 

0.5103 - 0.2858i 

0.5114 - 0.2854i 

0.5119 - 0.2863i 

2 

0 

0.7346 - 0.0949i 

0.7348 - 0.0949i 

0.7348 - 0.0949i 

2 

1 

0.7206 - 0.2870i 

0.7210 - 0.2869i 

0.7211 - 0.28711 

2 

2 

0.6960 - 0.4844i 

0.6953 - 0.4855i 

0.6892 - 0.4834i 

3 

0 

0.9343 - 0.0954i 

0.9344 - 0.0954i 

0.9344 - 0.0954i 

3 

1 

0.9233 - 0.2876i 

0.9235 - 0.2876i 

0.9235 - 0.2876i 

3 

2 

0.9031 - 0.4835i 

0.9026 - 0.4840i 

0.9026 - 0.4840i 

3 

3 

0.8759 - 0.6835i 

0.8733 - 0.6870i 

0.8733 - 0.6870i 

4 

0 

1.1315 - 0.0956i 

1.1315 - 0.0956i 

1.1315 - 0.0956i 

4 

1 

1.1224 - 0.2879i 

1.1225 - 0.2879i 

1.1225 - 0.28791 

4 

2 

1.1053 - 0.4828i 

1.1050 - 0.4831i 

1.1050 - 0.4831i 

4 

3 

1.0817 - 0.6812i 

1.0798 - 0.6830i 

1.0798 - 0.6830i 

4 

4 

1.0530 - 0.8828i 

1.0485 - 0.8891i 

1.0485 - 0.8891i 

5 

0 

1.3273 - 0.0958i 

1.3273 - 0.0958i 

1.3273 - 0.0958i 

5 

1 

1.3196 - 0.2881i 

1.3196 - 0.28811 

1.3196 - 0.28811 

5 

2 

1.3048 - 0.4824i 

1.3045 - 0.4826i 

1.3046 - 0.4826i 

5 

3 

1.2839 - 0.6795i 

1.2826 - 0.68051 

1.2826 - 0.6805i 

5 

4 

1.2582 - 0.8794i 

1.2548 - 0.8832i 

1.2548 - 0.8832i 

5 

5 

1.2284 - 1.0821i 

1.2221 - 1.09151 

1.2221 - 1.0915i 


The AIM procedure can now be applied to Eq. fH3l) and the result obtained after 150 
iterations are tabulated in Table HI One can see that the results are dehnitely better than 


the third order WKB and are comparable with the sixth order WKB and the prony results 
in Ref. jl^ . 
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V. ABSORPTION PROBABILITIES 


For the low energy regime analytic formula of the absorption probability in the scattering 


of gravitinos can be obtained using the Unruh approximation method 


291. On the other 


hand, for the more general energy regimes we shall use the WKB approximation, especially 
the one developed by Iyer and Will 141. 

To implement the Unruh method we begin with the potential in Eq. (jST]). We shall 
consider three regions: (I) Near horizon where /(r) 0, (II) Central region where the 

potential is much larger than the energy with U(r) ^ u, and (III) Far from the black hole 
where /(r) —)■ 1. Approximated solutions are obtained for these regions separately and 
then the unknown coefficients are determined by matching these solutions at the boundary 
regions. 

In the hrst region, /(r) —)■ 0, Eq. (I35|) becomes 

= (46) 

With the ingoing boundary condition near the event horizon, the solution goes like 


$1 = Aie~ 


(47) 


Next, in the region U(r) 3> cu, Eq. fl3^ becomes 


-^ + w 

dr^ 


A _ pp) = 0. 


Dehning H as 




Eq. 


transforms into a hrst-order differential equation 


The solution of Eq. fl50|) is given 

H = Bn 


d 

— + W\ H = Q. 
dr^ 


fj + h-yj 


yj 


j + \ + yj 


(48) 


(49) 


(50) 


(51) 


Inserting Eq. (15T]) into Eq. 
the solution goes like 

$11 = A 


, we will have a hrst order differential equation of and 


II 


i + yjy^^ fj + " 


yj 


yj 


j +1 + yj 


+ , 


(52) 
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where is 


$ = 


1 + v7 


i+5 


3 + \- vT 


. 1 ~ vT J Vi +1 + v7 

In the region where /(r) —)■ 1, Eq 


1 A-/7 


(ir^ 


$iii - 


/ Vi + yj 

becomes 

{f - J) 


2j+l 


i +1 + \/7 
i +1 - v7 


rfr' 


(53) 


— cu 


<hTTT = 0 . 


(54) 


Note that in this region r* ~ r, the solution will be expressed in terms of Bessel functions 
as 

*hiii = Alll^/rJj {ur) + Biii^/rNj (cnr) . (55) 

Letting the incoming part of $111 to have unit amplitude at r ^ cxo, we have the relation 
between Am and Bm which reads 


Am + iBm = V^StixJ . 


(56) 


Matching the solutions in regions I anc 
can obtain the absorption probability as 


II and also the solutions in regions II and III, one 


29 


30| 




(57) 


for small u, where 


C = 


3 + 


1 ; ’ 


(58) 


2'wr(j + i)y j 

and r is the gamma function. 

To evaluate the absorption probabilities for the whole energy range we adopt the WKB 
approximation. To implement the method it is convenient to change variable to a; = cur and 
to take Q {x) = — V^ where Eq. fl35l) then becomes 


dxl 


+ Q $1 = 0 . 


(59) 


For energy <^V the low energy absorption probabilities are given by the first order WKB 
approximation, which reads M| 


\Aj (cu) 1^ = exp 


a™/ ^_ 


(60) 
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where xi and X 2 are turning points where Q{xi,X 2 ) = 0 or i/ {xi,X 2 ) = uj"^, for a given 
energy u with potential V. 

When ~ fo the formula in Eq. flhOp will no longer be appropriate as the exponential goes 
to inhnity. In this energy regime a suitable method to derive the absorption probabilities is 
by using the third order WKB approximation of Iyer and Will [l^. Using the same notation 


as in Ref. 


where 


32| . the absorption probability can be expressed as 

lAi (uj) p = ^ 


1 g2S(a;) ’ 


(61) 


S{uj) = -zl + 


Oc! — —0, 


64 


16 


+7r/c^/^ 

+7r/c“^/^ 

-TTfc-'/' 


/1155 . 315^21. 35 ^2 35^ ^ 5 > g 

-be — + -zr^b^ + —^ 3^5 ~ 


V2048 


256 


128 


64 


32 


3. 7,2 

■Oa - bn 


16 
/1365, 


64 




525 .. 85 . 85,^2 95 ,^,^ 

-6064 - Ija -64 “ 1 “ —^3^5 

256 ^ 128 ^ 128 ^ 64 ® 


25, 


, (62) 


2048 ^ 256 128 ^ 128 ^ 64 32 

where Zq, bn, and k are defined by the components of the Taylor series in expanding Q{x) 
near its peak at xq, 


Q — Qo + ^ ~r 

nl 


71=3 


rQ 

dx'^ 


= k 


+ 


71=3 


(63) 


and 


z = X 


Xq 


k = -Q 
2 ^ 


// 

0 ) 


^2 ^ _9^0 

Q'o' 


brt. - 


n\Q 


d^Q 

dx'^ 


(64) 


The subscript 0 represents the maximum of Q and the primes denote derivatives. 

The absorption probabilities given by Eq. (16T1) from j = 3/2 to 11 / 2 , and are presented in 
Fig. |2j They are plotted as a function of energy u. As the parameter j is increased from 3/2 
to 11 / 2 , the absorption probabilities shift to the right, and all the absorption probabilities 
approach 1 as the energy becomes large, that is, in the high energy regime, as expected. 
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FIG. 2: Absorption probabilities by the WKB approximation for j = 3/2 to 11/2. 

VI. CONCLUSIONS AND DISCUSSIONS 


In this paper we have introduced a novel method to deal with the spin-3/2, or the grav- 
itino, helds in spherically symmetric spacetimes. Our method can be readily extended to 
higher dimensions, in contrast to the Newman-Penrose formalism which is more specihc to 
four dimensions. In this method we considered the eigenspinor-vectors on higher-dimensional 
spheres. Actually the eigenspinor-vectors of the Dirac operator in higher-dimensional max¬ 
imally symmetric spaces is an interes ting topic of its own. Along the same lines as the 
work done by Camporesi and Higuchi 33|], for the case of the eigenspinors, one can obtain 
the spinor-vector eigenfunctions, their eigenvalues and degeneracies on spheres as well as on 
hyperbolic spaces. Work in this direction is in progress. 

With the help of these eigenspinor-vectors it is possible to consider spherically symmetric 
spacetimes in higher dimensions. As we have mentioned earlier, with the addition of the 
spin-3/2 helds, we can establish the QNM frequencies and the greybody factors in Hawking 
radiation processes for all spins in higher dimensional spherically symmetric black holes. 
This will be helpful in understanding the production and the evolution of possible collider 
produced TeV black holes |^. 

Other than the higher dimensional Schwarzschild black holes one could apply our method 
to other spherically symmetric spacetimes like the Reissner-Nordstrom as well as the AdS- 
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and dS-black holes. For the charged black hole cases, especially the extremal ones, we could 
examine the behaviors of supersymmetric black holes. For example, in 0 it is found that 
the asymptotic QNM frequencies are the same for spin-0, 1/2, 1, 3/2 and 2 helds in the 
extremal case in four dimensions. It is therefore interesting to see if this holds for extremal 
black holes in higher dimensions. 

Our method should also be useful in the study of properties of AdS black holes. These 
will be important in relation to the AdS/CFT correspondence 3^ to understand various 
hnite temperature behaviors of the boundary conformal held theory 3^. We can examine 


the gravitino perturbations for the whole space outside the event horizon 


36|, not just for 


the near horizon region. We shall also address this problem in a future work. 
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Appendix A: Eigenspinors and eigenspinor-vectors on the 2-sphere 

The metric of a 2-sphere can be taken as 

= sin^ 6'2 dOl + dO^ , (Al) 

where the zweibein are dehned as 

e/ = diag(sin 6 ' 2 ,1) ; e/ = diag(l/sin 6 * 2 ,1) . (A 2 ) 

Note that we have omitted the overbars for simplicity in this appendix. 

The Dirac operator is given by 

7“V,^ = 7V(5/. + r^)^- (A3) 
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In this metric the non-vanishing Christoffel symbols and spin connections are 


= cot 62 ; Fg" = - sin 6*2 cos 6^2 , 


e »2 


0102 


0101 


and 


(A4) 


F 01 = - cos 62 


F^a — 0 . 


(AS) 


The Dirac matrices are chosen to be 




a 


7 ®^ = 


, sin 62 , 

where cr* are the Panli matrices. 

Snppose we write the Dirac operator eigenspinor equation as 


(A 6 ) 




-o-^^^ 50 i + ^02 + 7 cot 02 

sm O 2 V 2 


ipx = iXi>x , 


(A7) 


where iX is the eigenvalue. Consider first the eigenvalue equation for dg^^, 


dg^Xm = imxrn , 


(A 8 ) 


that is, Xm ~ Note that for spinors one should get a sign change for a 27r rotation in 

9i. Therefore, the eigenvalues of m should be half-integers. 


1 3 5 

m = ±-, ±-, ±-, 
2 ’ 2 ’ 2 ’ 


Going back to the eigenvalue equation for we write 


(A9) 


i^x = 


Ax{O2)Xm{01 


Bx{e2)xm{e^) 

Putting this into the eigenvalue equation, we have the following set of equations. 


(AlO) 


m 


dg^ + 7 cot 92] Ax + ^-^Ax = iXBx , 

2 J sm 6*2 

dg^ + - cot 92 ] Bx- -t—^Bx = iXAx . 

2 J sm 02 

These can be turned into second order equations. For example, for Ax we have. 


(All) 


<902 + X cot 02 -: 


m 


sin 02 


^02 + o cot 02 + w 


m 


sin 02 


Ai = -AM> 


(A 12 ) 
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This equation can be further simplified by dehning 


4a = (sin « 2 ) “A ■ 


(A13) 


Eq. flA12D then becomes 

(fux 

del 

d^Uy 

with the solution |3 


+ A 


+ 


2 rn? m cos 6*2 A 

^ 2 T 1 — 0 

sm 62 sm 62 J 

1 _ 4 (^ + 1)2 ^ 1-4 (^-1)2 ^ 


16 sin 2 ^ 


16 cos^ Y 


ma = 0 , 


U\ = sm 


02 A f 


( fj \ 2 

COS^j Pi“’^)(cos02) 


Pi°‘’^\cos 62 ) are the Jacobi polynomials with 


a = 


m 


(9 = 


m — 


; n = ±A — 


CK T /9 T 1 


n = 0,1, 2,... for to be a polynomial, so that it is convergent in 0 < 02 < tt. Hence, 

for the values of m given in Eq. flAOp . the eigenvalues of the eigenspinors are 


(A14) 


(A15) 


(A16) 


A = ±1,±2,±3, 


(A17) 


Next we consider the eigenspinor-vectors of the Dirac operator. 


(A18) 

on the 2-sphere. From the eigenspinor 'i/jx, one can construct the two independent 
eigenspinor-vectors on using the linear combinations of and 

= V^V'A + ai7/.'0A ; = V^-^A + a27M'*/'A , (A19) 

where 

ai='- (-A + \/A2 - 1 ) ; ei = \/A2-l , (A20) 

02 = ^ (-A - \/A2 - 1 ) ; 5 = -VA2 - 1 . (A21) 

They are orthogonal to each other in such a way that 

/" dD = 0 . (A22) 
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fl 2) 

Note that they do not obey the condition 7^'^^ ’ — 0. Indeed, 

7^(V^V'A + ai,27M^A) = (*A + 2ai,2)^/'A 

= ±WX^ - 1 ipx ■ (A23) 

These vanish only for the zero modes, ^1,2 = 0, which we will exclnde. Also, we have 

= -^\/A2 - 1 (±A + \/A2 - 1) V'A , (A24) 

which are also non-vanishing. 

Hence, from the eigenvalues of the eigenspinors in Eq. flA17D . we have the eigenvalues of 
the eigenspinor-vectors on the 2-sphere, 

6,2 = ±\/A2 - 1 = ±a/3, ±x/8, ±x/15, ..., (A25) 

(1 2 ) 

while the two eigenmodes V'm = V^V'a + Oigd/^t/'A are linear combinations of V^V'a and 
7/.'0 a- 
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